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Abstract 

We simulate several models of random curves in the half plane and numerically com- 
pute the stochastic driving processes that produce the curves through the Loewner equa- 
tion. Our models include models whose scaling limit is the Schramm-Loewner evolution 
(SLE) and models for which it is not. We study several tests of whether the driving 
process is Brownian motion, as it is for SLE. We find that testing only the normality of 
the process at a fixed time is not effective at determining if the random curves are an 
SLE. Tests that involve the independence of the increments of Brownian motion are much 
more effective. We also study the zipper algorithm for numerically computing the driving 
function of a simple curve. We give an implementation of this algorithm which runs in 
a time 0(N 1,S5 ) rather than the usual 0(N 2 ), where N is the number of points on the 
curve. 
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1 Introduction 



The Loewner equation provides a means for encoding curves in the upper half plane that do not 
intersect themselves by a real- valued function. Let j(t) be such a simple curve with < t < oo. 
Let H denote the upper half of the complex plane, and let 7[0, t] denote the image of 7 up to 
time t. Then H \ 7[0,t] is a simply connected domain. So there is a conformal map g t from 
this domain to H. If the curve is suitably parametrized and gt is suitably normalized, then g t 
satisfies the differential equation 



dgt(z) _ 2 
dt g t {z)-U t 



9o{z) = z (1) 



for some real valued function U t on [0, 00). The function U t is often called the driving function. 

If our simple curve in the half plane is random, then the driving function U t is a stochastic 
process. Schramm discovered that if the scaling limit of a two-dimensional model is conformally 
invariant and satisfies a certain Markov property, then this stochastic driving process must be 
a Brownian motion with mean zero [H]. The only thing that is not determined is the variance. 
Schramm named this process stochastic Loewner evolution or SLE; it is now often referred to 
as Schramm-Loewner evolution. 

Many critical two-dimensional models from statistical mechanics and probability satisfy 
these properties or are believed to satisfy them, and so should be SLE for some k (the parameter 
that determines the variance). These include the loop-erased random walk [TTJ [19], the self- 
avoiding walk [10], interfaces in the critical Ising model [16] . the Gaussian free field [17], critical 
percolation [13 [5], and uniform spanning trees [UJ. More recent work has considered whether 
other models have random curves that are described by SLE. The possibility that domain walls 
in spin glass ground states are SLE curves was studied numerically both by Amoruso, Hartman, 
Hastings, and Moore [JJ and by Bernard, Le Doussal, and Middleton [4j. Bernard, Boffetta, 
Celani and Falkovich considered simulations of certain isolines in two-dimensional turbulence 
[2] and surface quasi-geostrophic turbulence [3]. 

In this paper we consider models which are definitely not SLE. They are based on well 
known lattice models - the loop-erased random walk (LERW), the self-avoiding walk (SAW) 
and the critical percolation exploration process. We distort these models by shrinking the 
random curves slightly in the vertical direction but not in the horizontal direction. In other 
words we apply a non-conformal transformation to the curves. Without distortion these models 
are all proven or conjectured to have a scaling limit given by SLE. 

One way to test if a model of random curves is SLE is to compute its stochastic driving 
process and see if it is Brownian motion. In this paper we simulate these distorted models, 
numerically compute their stochastic driving process, and then test if they are Brownian mo- 
tions. We also do this for the models without distortion. Our goal is to see how well one can 
determine whether or not a model is SLE by studying this stochastic driving process and to 
compare various methods for testing if the stochastic driving process is a Brownian motion. 

Another goal of this paper is to study the algorithm for computing the driving function of 
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a given curve. The standard implementation of the "zipper algorithm" for doing this requires 
a time 0(N 2 ) where N is the number of points on the curve. We present an implementation 
that runs in a time 0(N P ) with p approximately 1.35. This implementation uses the same idea 
used in [6j to simulate SLE curves quickly. We also study the difference in the driving function 
found using "tilted slits" versus "vertical slits" and the effect of the number of points used on 
the curve to compute the driving function. 

2 Distorted models 

We study the stochastic driving function of three models. We refer to them as distorted models. 
For A > we define a non-conformal map on the upper half plane by <p(x, y) = (x, Xy). Given a 
model that produces random curves 7 in the upper half plane, we consider the random curves 
</> o 7. In other words, we stretch the curve by a factor of A in the vertical direction, but do not 
stretch it in the horizontal direction. 

We apply this distortion to the loop-erased random walk, the self-avoiding walk and perco- 
lation interfaces. For all three models we consider the chordal version of the model in which 
the random curve lies in the upper half plane and goes from the origin to 00. For the LERW 
this means we take a half plane excursion and loop erase it. For the SAW this means we use 
the uniform probability measure on nearest neighbor walks with a fixed number of steps which 
begin at the origin and lie in the upper half plane. For percolation we consider site percola- 
tion on the triangular lattice in the upper half plane with boundary conditions which force an 
interface to start at the origin. Details of the definitions of these models and the parameters 
used in the simulations may be found in the appendix. 

There is no simple relation between the driving function for the curve 7 and the driving 
function for the distorted curve <fi o 7. We study the driving function for the distorted model 
as follows. We generate samples of the LERW, SAW or percolation interface and then apply 
the distortion map 0. Then we compute the driving function of the distorted curve. The result 
is a collection of samples of the stochastic driving process of the distorted model. We then do 
various statistical tests to see if this process is a Brownian motion. 

We denote the driving process by U t . All the models are invariant under reflections about 
the vertical axis. Hence E[U t ] = 0. We begin by plotting the variance E\Uf] as a function of 
t. We should emphasize that in the scaling limit, all the models have a scaling property which 
implies that E[U 2 ] is a linear function of t, even if the scaling limit is not an SLE. So this 
does not test whether U t is a Brownian motion. It only provides an estimate of k where k is 
the slope of the function t — > E[Uf]. 

The first statistical test is to see if the distribution of an individual U t is normal. We use the 
Kolmogorov-Smirnov test. This test is based on the fact that for a continuous random variable 
Y, if F is the cumulative distribution of Y, then F(Y) is uniformly distributed on [0, 1]. Let 
Yi,Y 2 , ■ ■ ■ ,Y N be N observations of the random variable Y, and let K(i) < K( 2 ) < • • • Y(jv) be 
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these numbers arranged in increasing order. Then the statistic is 

k + - 1 

D= max \F(Y (k) ) -4 1 + — (2) 

i<k<N l y ( >J N 1 2N y ' 

(This formula may appear different from that found in most statistics texts, but it is in fact 
the same.) Under the null hypothesis that Y has the distribution F, the limiting distribution 
of \fND as N — > oo is known. For example, P(\/ND > 1.36) is approximately 5%. So if we 
compute this statistic for an individual U t with F equal to the cumulative distribution for a 
normal random variable with mean zero and variance ut, and find that the value of y/~ND is 
large (say larger than 1.36), then we conclude that the distribution of U t is not this normal 
distribution. We perform this Kolmogorov-Smirnov test for two values of the time, T and T/2. 
The value of T as well as other parameters used in the simulations are given in the appendix. 
The results are shown in the various tables in the columns labeled D(T/2) and D(T). 
Our next test involves the independence of increments of Brownian motion. Let 

X 1 = U T /2, X 2 = U T — U T /2 (3) 

If U t is a Brownian motion, then X 1 and X 2 are independent and have mean zero. So their 
product XiX 2 has mean zero. We test the hypothesis that its mean is zero. If U t is a Brownian 
motion, then the variance of X\X 2 can be found. Denote it by a 2 . The statistic we use, Z, 
is simply the sample mean for X\X 2 divided by a/y/N. If N is large and Ut is a Brownian 
motion, then the distribution of Z is close to that of the standard normal. This test is shown 
in the tables in the column labeled Z. 

For our last three tests we let < ti < t 2 < ■ ■ ■ < t n = T. We use equally spaced U in 
our tests, but one could consider non-uniform choices. If U t is a Brownian motion, then the 
increments 

X j = U tj - U tj _, (4) 

are independent and each is normal with mean zero and variance n(tj — tj-i). We test this 
joint distribution with a x 2 goodness-of-fit test. The idea is to divide the possible values of 
(X 1: X 2 , ■ ■ ■ , X n ) into m cells and count the number of samples that fall into each cell. Under 
the hypothesis that the Xj are independent and normally distributed, we can compute the 
expected number of samples that fall into each cell. Let Oj be the number of samples in cell j, 
and Ej the expected value of this random variable under the hypothesis that Ut is Brownian 
motion. The test is then based on the statistic: 

x2 = £(o l _E£ (5) 



J'=l 
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Under the hypothesis that U t is a Brownian motion, if m is large the distribution of this statistic 
is approximately the \ 2 distribution with m — 1 degrees of freedom. 

We consider three choices of the cells. In the first choice we let n — 10, and use only the signs 
of the Xj to define the cells. Thus there are 2 10 = 1024 cells. We denote the statistic in this case 
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Figure 1: The points are a histogram for the density of Ut/VT for the LERW, SAW 
and percolation with distortion A = 0.95. The curves are the density of a normal 
distribution with variance k where k is determined from the least squares fit. 

by xi- The second choice of cells is similar - we use only the signs of the increments but with 
n — 5. So there are 2 5 = 32 cells. This statistic is denoted by x\- The third choice uses just two 
increments. For each of the increments X\ and X2 we look at which quartile it falls in. We let 
q be the constant such that under the hypothesis that U t is Brownian motion, P(Xi > g)=l/4, 
and then divide the possible values of each into the intervals (—00, — q], [—q, 0], [0, q), [q, 00). 
Thus there are 16 cells. This final statistic is denoted by xt- Note that the statistics x\ an d 
x\ have the advantage that they do not involve the value of k. For x\ we need a value for k to 
compute q. 

We study three values of the distortion parameter, A = 0.9,0.95 and 1, for each of the 
three lattice models. Thus there are a total of nine models considered. (We have run the same 
simulations for A = 1.05 and 1.1, but do not present these results. They are consistent with 
the results that we do present.) We generate 10 5 samples for each of the nine cases, compute 
the driving function U t for each sample, and record the value of U t at ten equally spaced times, 
t = T/10, 2T/10, • • ■ , T. We compute the variance of U t at each of the ten times and then do 
a least squares fit to estimate k, the slope of t — > E[Uf]. The results are shown in table [TJ 

For each of the three models we present the results of our statistical tests in a table. 
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MdHpI 


\ 


K 


±J±1/±VV V 


n Qfi 


O ooi r _i_ n m r\e 
Z, .OO LO _T_ U.U1UU 


LERW 


0.95 


2 1709 ± 0094 


LERW 


1.00 


2.0008 ± 0.0093 


SAW 


0.90 


3.0744 ± 0.0138 


SAW 


0.95 


2.8414 ± 0.0108 


SAW 


1.00 


2.6686 ± 0.0132 


percolation 


0.90 


6.9446 ± 0.0287 


percolation 


0.95 


6.4422 ± 0.0311 


percolation 


1.00 


6.0404 ± 0.0265 



Table 1: The estimate of k using a weighted least squares fit for E[U^} as a function 
of t. The error bars are two standard deviations. 



We perform the statistical tests for the first iV samples of the 100K samples, using N = 
5K, 10K, 20K, 50K and 100K. In the tables we do not give the value of the statistic, but 
rather the corresponding "p-value." The p-value is defined as follows. Consider a one-sided 
statistic such as the Kolmogorov-Smirnov D statistic. Suppose that in our test the value of 
the statistic is x. Then the p-value is the probability P(D > x) under the null hypothesis. For 
a two-sided statistic such as Z the definition is modified in the obvious way. A small p-value 
(less than a few percent) indicates that the value of the corresponding statistic is very unlikely 
under the hypothesis that Ut is Brownian motion, and so we should reject the hypothesis that 
Ut is Brownian motion. 

These nine tables contain a lot of numbers, but they clearly show the following. For A = 1 
the p-values in the table for all three lattice models are not small, and so our tests do not 
indicate that we should reject the hypothesis that U t is a Brownian motion. This is as it should 
be. When A = 1 the models should have a scaling limit given by SLE, and so U t should indeed 
be Brownian motion. 

For the models that are not SLE, A = 0.9, 0.95, we first consider the two Kolmogorov- 
Smirnov tests that U t is normal. The corresponding p-values are not typically small, and this 
test is completely ineffective at indicating that U t is not Brownian motion, even when we use 
10 5 samples. In figure [T] the data points are histograms for the density function of Ut / VT for 
the LERW, SAW and percolation with A = 0.95. The curves are the density functions for the 
normal distribution with variance k where k is determined from our least squares fit. As the 
Kolmogorov-Smirnov test showed, the data points are fit very well by the normal curves. (For 
A ^ 1, it is easy to show that the distorted model is not SLE, and so Ut is not a Brownian 
motion. However, this does not rule out the possibility that the Ut are normal even for the 
distorted models.) 

The other four tests involve the independence of the increments. The test based on the 
mean of the product of two independent increments sometimes indicates correctly that U t is 
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A 


N 


D(T/2) 


D(T) 


z 


y 2 

An 


Ab 


y 2 

Ac 




5 000 

W .WWW 


018695 

W . W 1UU JtJ 


223875 

W . — — Ww 1 W 


012146 


306555 


6541 53 


139523 




10 000 

1U .www 


605601 


562645 

W. W W i 5_/ 


012238 


000614 

W.WWWW X t: 


000666 

w. wwwwww 


352671 


0.90 


20,000 


0.578409 


0.860558 


0.029349 


0.000000 


0.000041 


0.820988 




50,000 


0.213470 


0.308500 


0.078267 


0.000000 


0.000000 


0.000120 




100,000 


0.222762 


0.253184 


0.003394 


0.000000 


0.000000 


0.000000 




5,000 


689088 


567792 


021239 


913735 


598610 

W. (JJUU X w 


766046 

W. 1 W W W i w 




10 000 

1U ■ w w w 


562020 

W .WW^jW^jW 


898586 

W. UJUl/UU 


055965 

W. V 7 :.J 'J •../ V.) "J 


423453 


350098 

W. Uuuu JU 


593728 


0.95 


20,000 


0.648585 


0.638202 


0.064346 


0.793825 


0.318862 


0.865687 




50,000 


0.486077 


0.178479 


0.010384 


0.078510 


0.000173 


0.850495 




100,000 


0.234619 


0.004251 


0.000812 


0.000080 


0.000000 


0.501958 




5,000 


0.322835 


0.909183 


0.607075 


0.599810 


0.281570 


0.332253 




10,000 


0.902220 


0.499358 


0.274931 


0.118588 


0.593618 


0.820108 


1.00 


20,000 


0.763553 


0.856747 


0.240621 


0.167155 


0.227972 


0.322061 




50,000 


0.840997 


0.746111 


0.376101 


0.695880 


0.735182 


0.702162 




100,000 


0.949877 


0.934915 


0.448990 


0.305246 


0.257668 


0.509009 



Table 2: The distorted LERW. A is the amount of distortion, with A = 1 being no 
distortion. iV is the number of samples used. The other six columns give the p-value 
of six different statistics used to test if the driving process is Brownian motion. See 
the text for details. 



not Brownian motion, but it is not very powerful. By contrast y 2 and y 2 are quite effective 
at correctly indicating when U t is not a Brownian motion. For A = 0.9 these tests correctly 
indicate Ut is not a Brownian motion with sample sizes on the order of ten to twenty thousand. 
For A = 0.95 these tests need on the order of a hundred thousand samples, but they are the 
only tests to correctly indicate that Ut is not a Brownian motion for this amount of distortion. 
The final statistic y 2 is only sometimes effective. 

3 A faster zipper 

We briefly describe the standard method for computing the driving function of a simple curve 
7. Let g s be the conformal map which takes the half plane minus 7[0, s] onto the half plane, 
normalized so that for large z 

2t 1 

g s ( Z )=Z + - + 0(^), (6) 

The coefficient 2t depends on s and is the half-plane capacity of 7[0,s]. The value of the 
driving function at t is U t = g s ( r y(s)). Thus computing the driving function essentially reduces 
to computing this uniformizing conformal map. We will describe the "zipper algorithm" for 
doing this [5JQ2]- Another approach to computing the driving function may be found in |18j . 
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A 


N 


D(T/2) 


D(T) 


z 


y 2 

Aa 


Ab 


y 2 

Ac 




5 000 


076330 


D 174654 


5381 22 


218454 


170313 


903646 




10 000 


471555 


343790 


824175 


100552 


019827 


803437 


0.90 


20,000 


0.246740 


0.321011 


0.905811 


0.003206 


0.000008 


0.464170 




50,000 


0.165784 


0.312940 


0.013581 


0.000000 


0.000000 


0.004110 




100,000 


0.251134 


0.200176 


0.000124 


0.000000 


0.000000 


0.000000 




5,000 


408348 


350261 


933527 


027830 


037820 


883070 




10,000 


136263 


888704 

U. UUU 1 V/ i 


681952 


775130 


704539 


931479 


0.95 


20,000 


0.805490 


0.663008 


0.983501 


0.682962 


0.797951 


0.104660 




50,000 


0.483299 


0.323554 


0.882502 


0.162307 


0.030086 


0.959757 




100,000 


0.313017 


0.158268 


0.794145 


0.000018 


0.009703 


0.534966 




5,000 


0.456948 


0.371554 


0.636449 


0.618423 


0.126116 


0.424437 




10,000 


0.557806 


0.267715 


0.821561 


0.200096 


0.145614 


0.609878 


1.00 


20,000 


0.933089 


0.552806 


0.514615 


0.793825 


0.253454 


0.370379 




50,000 


0.956776 


0.852597 


0.588566 


0.685205 


0.275748 


0.161052 




100,000 


0.501460 


0.870474 


0.219944 


0.028353 


0.079678 


0.579479 



Table 3: The p-values of the distorted SAW. 



We find it more convenient to work with the conformal map 

h a (z)=g 8 (z)-U 9 (7) 

It maps H\7[0, s] onto EI and sends the tip 7(5) to the origin. The value of the driving function 
at s is minus the constant term in the Laurent expansion of h s about 00. From now on we 
work with this normalization for our conformal maps. 

Let z , zi, ■ ■ ■ , z n be points along the curve with z = 0. In our applications these are lattice 
sites. The zipper algorithm finds a sequence of conformal maps hi, i = 1,2, • • • ,n such that 
hk o hk-i o ■ ■ ■ o hi approximates the conformal map for the curve up to site z^. Suppose that 
the conformal maps h\, hi, • • ■ , have been defined so that h^ o /i fc _ 1 o • • • o hi sends HI \ 7 to 
HI where 7 is some curve that passes through z , Zi, • • • z^. In particular z k is mapped to the 
origin. Let 

w k+ i = h k o h k -i o ■ ■ ■ o hi(z k+1 ) (8) 

Then w k+ i is close to the origin. We define h k +i to be a conformal map with the appropriate 
normalizations that sends HI \ 7 fc+1 to HI where ^ k+ i is a short simple curve that ends at w k+1 . 
The key idea is to choose this curve so that h k+ i is explicitly known. The two choices we will 
use are "tilted slits" and "vertical slits." 

Let 2Atj be the capacity of the map hi, and AUi the final value of the driving function for 
hi. So 

h i (z) = z-AU i + ^ + 0(\) (9) 



8 



A 


N 


D(T/2) 


D(T) 


z 


Y 2 
Aa 




y 2 

Ac 




5 000 

W .WWW 


26881 1 


707407 


035275 


463308 


n 440588 

W . tit: UuUU 


068530 

w . w WO WOW 




10 000 

1U .WWW 


568890 


306491 


043429 

W • V 7 i W i Z_ .../ 


387697 


125325 


286280 


0.90 


20,000 


0.113505 


0.640373 


0.192064 


0.002475 


0.000123 


0.387076 




50,000 


0.073017 


0.749263 


0.000274 


0.000000 


0.000000 


0.007620 




100,000 


0.000694 


0.015332 


0.000005 


0.000000 


0.000000 


0.000000 




5,000 


667251 

W . W W l i-i\J X 


576948 

W . W l W t7 


881721 

W • O W -L 1 ill 


978563 


820268 

w. u^u^uu 


261246 




10,000 


831016 

W.tJWXWXW 


921401 


964028 


586700 

\J • ^JkJKJ I WW 


601949 

W.WW IJ^J 


198512 


0.95 


20,000 


0.827859 


0.938272 


0.703387 


0.580099 


0.123268 


0.124509 




50,000 


0.504299 


0.913135 


0.465606 


0.033206 


0.000071 


0.291786 




100,000 


0.084361 


0.515339 


0.259067 


0.000000 


0.000000 


0.043819 




5,000 


0.075710 


0.685162 


0.242412 


0.884633 


0.296115 


0.657406 




10,000 


0.882136 


0.304884 


0.618193 


0.267297 


0.246049 


0.766046 


1.00 


20,000 


0.271328 


0.037469 


0.133359 


0.088622 


0.587819 


0.060798 




50,000 


0.100113 


0.008876 


0.184139 


0.256426 


0.996808 


0.283681 




100,000 


0.333954 


0.059206 


0.518117 


0.287868 


0.779478 


0.149034 



Table 4: The p-values of the distorted percolation interface. 



Then 

2t 1 

h k o /ifc_i o • • • o ^(z) = 2 - C/ t + — + 0{— ) (10) 

where 

k 

i=i 

f/ t = ^Af/, (12) 

i=l 

Thus the driving function of the curve is obtained by "adding up" the driving functions of the 
elementary conformal maps hi. 

We now consider the two particular types of maps we use for hk+i- For tilted slits, 7^+1 is 
the line segment from the origin to Wk+i- There is no explicit formula for hk+i in this case, but 
there is a formula for its inverse: 

h£ 1 (z) = (z + x l ) 1 - a (z-x r ) a (13) 

where xi, x r > 0. It maps the half plane onto the half plane minus a line segment which starts at 
the origin and forms an angle a with the positive real axis. The interval [—2^,2;,.] gets mapped 
onto the slit. We must choose x\ and x r so that hk+i satisfies our normalization conditions. 
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In particular, /i fc _f 1 must send the origin to the tip of the line segment, i.e., Wk+i- Tedious but 
straightforward calculation shows if we let w k +\ = r exp(ict7r) , then 



1 — a \ a 



Xl = r {—) ' Xr = r {l^) (M) 
The changes in the driving function are given by 

At = ^r 2 a l - 2a (l-a) 2a -\ AU = r(l - 2a)a~ a (l - a)- {1 - a) (15) 

For vertical slits we take 7^+1 to be the vertical line segment from the real axis to Wk+i- 
Note that since this 7^+1 does not start at the origin, this method does not approximate the 
original simple curve 7 by another simple curve. Instead, the domain of the conformal map 
that we construct to approximate h s is of the form H minus a set that is more complicated than 
a simple curve. This may make the reader nervous, but we will see in the next section that we 
get essentially the same driving functions using tilted slits or vertical slits. The conformal map 
that removes this vertical line with our normalizations is 

h k+ i(z) = iy/-(z - x) 2 - y 2 (16) 

where w^+i = x + iy and the branch cut for the square root is the negative real axis. The 
changes in the driving function are given by 

At = ^y 2 , AU — x (17) 

A comment on terminology is in order. We use "zipper algorithm" to refer to all the various 
algorithms we can get from different choices of the curve 7^+1 • Marshall and Rohde [13] use 
"zipper" to refer only to the choice using tilted slits. 

The number of operations needed to compute a single w k +i is proportional to k. So to 
compute all the points w k +i requires a time 0(N 2 ). To do better we must avoid evaluating the 
fc-fold composition in (jSj) every time we compute a Wk+i- We begin by grouping the functions 
in (jSJ) into blocks. We denote the number of functions in a block by b. Let 

Hj = h jb o h jb -i o • • • o h(j-i)b+2 ° ^0-i)b+i (18) 

If we write k as k = mb + r with < r < b, then (jHJ) becomes 

Wk+i = h mb+r o h mb+r _i o • • • o h mb+ i o H m o H m -i o • • • o Hi{z k+ i) (19) 

Typically, the number of compositions in ( fl9i) is smaller than the number in (jHJ) by roughly a 
factor of b. The hi are relatively simple, but the composition Hj cannot be explicitly computed. 
Our strategy is to approximate the hi by functions whose compositions can be explicitly com- 
puted to give an explicit approximation to Hj. This allows us to compute the compositions in 
(ITS]) just once rather than every time we compute a w k . 
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Recall that hi is normalized so that hi(oo) = oo and h'^oo) = 1. It maps EI minus a simple 
curve which starts at the origin to H. Let h denote such a conformal map. Let r be the largest 
distance from the origin to a point on the curve. Then h is analytic on {z G EI : \z\ > r}. Note 
that h is real valued on the real axis. By the Schwarz reflection principle it may be analytically 
continued to {z G C : \z\ > r}. Moreover, it does not vanish on this domain. So if we let 
f(z) = h(l/z), then / is analytic in {z G C : \z\ < 1/r} and /(0) = 0, /'(0) = 1. The Laurent 
series of h about oo is just the power series of / about 0. For large z, h(z) is well approximated 
by a finite number of terms in this Laurent series. It will prove more convenient to work with 
a different series. 

Define h(z) = l/h(l/z). Since h(z) does not vanish on {z G C : \z\ > r}, h(z) is analytic 
in {z G C : \z\ < 1/r}. Our assumptions on h imply that h(0) = and h (0) = 1. So h has a 
power series of the form 

oo 

h{z) = ^ (20) 

with oi = 1. The radius of convergence of this power series is 1/r. Note that the coefficients of 
this power series are the coefficients of the Laurent series of 1/h. 

The primary advantage of working with the power series of h is its behavior with respect 
to composition. It is trivial to check that 

{hi o h 2 y = hi o ti 2 (21) 

Our approximation for hi(z) is to replace hi(z) by the truncation of its power series at order n. 
So 

1 



hi(z) 



hiil/z) 



J2 a J z ~ 

.i=i 



(22) 



For each hi we compute the power series of hi to order n. We use them and ( !21j) to compute 
the power series of Hj to order n. Let 1 / Rj be the radius of convergence for the power series of 
Hj. Now consider equation ( TL9l) . If z is large compared to Rj, then Hj(z) is well approximated 
using the power series of Hj. We introduce a parameter L > 1 and use this series to compute 
Hj(z) whenever \z\ > LRj. When \z\ < LRj, we just use (fl8j) to compute Hj(z). The argument 
of Hj is the result of applying the previous conformal maps to some Zk+i, and so is random. Thus 
whether or not we can approximate a particular Hj by its series depends on the randomness 
and on which Wk+i we are computing. 

We need to compute Rj. Consider the images of %-i)6, z^j-i^+i, ■ ■ ■ Zjb-i under the map 
Hj-ioHj-2°- ■ -oH\. The domain of the conformal map Hj is the half-plane EI minus some curve 
Yj which passes through the images of these points. The radius Rj is the maximal distance from 
the origin to a point on Yj. This distance should be very close to or even equal to the maximum 
distance from the origin to images of z^-iu, %_i)&+i, • ■ ■ under Hj_i o Hj_ 2 o • • • o Hi. So 
we take Rj to be the maximum of these distances. 
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Table 5: The time (in seconds) needed to unzip a SAW with N steps without using 
the power series approximation is shown in the second column (time 1) The time 
using the power series approximation is shown in the third column (time 2). The 
fourth column (factor) is the ratio of these two times. The block length used is in the 
last column. 



Our algorithm depends on three parameters. The integer b is the number of functions in a 
block. The integer n is the order at which we truncate the power series of the Hj. The real 
number L > 1 determines when we use the power series approximation for the block function. 

The improvement in the speed of the zipper algorithm from using our power series approx- 
imation is shown in table [5] and figure EJ In these timing tests we use a single SAW with 
one million steps. We time how long it takes to unzip the first N steps with and without the 
power series approximation. We do the computations using the power series approximation 
for different choices for the block length, namely b = 20, 30, 40, 50, 75, 100, 200, 300, and report 
the fastest time. The last column in the table indicates the block length that achieves the 
fastest time. As a rule of thumb, a good choice for the block length (at least for the SAW) is 
b = yN /4. The next to last column in the table gives the factor by which the use of the power 
series approximation reduces the time needed for the computation. These timing tests were 
done on a PC with a 3.4 GHz Pentium 4 processor. 

Without the power series approximation the time is 0(N 2 ). This is seen clearly in the log- 
log plot in figure [2] where the data for unzipping without the power series approximation is fit 
quite well by a line with slope 2. The data for unzipping using the power series approximation 
is fit by a line with slope 1.35. This indicates that the time required when the power series are 
used is approximately 0(iV L35 ). 
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Figure 2: The points are the time (in seconds) needed to unzip a SAW with N steps 
with and without the power series approximation. The lines have slopes 2 and 1.35. 

4 Comparisons of computational methods 

Given a simple curve, there are several choices when we compute its Loewner driving function. 
For the explicit conformal map in the zipper algorithm we could use either the map that removes 
a vertical slit from the half plane or the map that removes a tilted slit. We could use the power 
series approximation to speed up the calculation as explained in the last section. We have a 
choice of how many points we use along the curve. In this section we study the effect of these 
various choices. 

For our study we generate a collection of forty self-avoiding walks with one million steps. 
For each walk we compute its driving function by three different methods. The first uses the 
vertical slit map and the power series approximation. The second uses the tilted slit map and 
the power series approximation. The third uses the vertical slit map without the power series 
approximation. To study the effect of using different numbers of points along the curve we 
do these computations on subsets of the million points on the SAW. The subsets are obtained 
by taking every nth point along the walk. For the vertical slit map with the power series 
approximation (the fastest of the three methods) we use n = 1,2, 5, 10, 20, 50, 100. So the 
curves being unzipped have from 1, 000, 000 to 10, 000 points on them. For the tilted slit map 
with the power series approximation we use n — 2, 5, 10, 20, 50, 100. (This method is slower 
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than that using vertical slits because of the need to use Newton's method in the computation of 
the tilted slit conformal map.) For the vertical slit map without the power series approximation 
(by far the slowest of the three methods) we use n — 5, 10, 20, 50, 100. 

Our choice of how to compare the driving functions computed by different methods for the 
same curve merits some discussion. One might measure the difference between two driving 
functions by computing the supremum norm or the L l norm of their difference over a bounded 
interval. We do not do this. The driving functions we are computing are approximations to 
Brownian motion sample paths. In particular, their slopes can be quite large. If you translate 
such a function by a small amount, the difference between the translated function and the 
original function can have a supremum norm or L 1 norm that is rather large. (Of course, the 
driving functions are continuous so these norms of the difference go to zero, but not linearly with 
the size of the translation.) Thus a small error in computing the capacity produces a relatively 
large error in these norms. Instead of using these norms, we compare driving functions by using 
only the last point on the driving function. We denote this last point by (T, Uj) ■ The capacity 
of the SAW is 2T and Ut is the image of the end of the SAW under the conformal map that 
uniformizes the half plane minus the SAW. 

We do not know the exact driving function of the SAW, so we treat the result of our 
computation using all one million steps (i.e., n — 1) with the vertical slit map and the power 
series approximation as the exact answer. We then compute relative errors for T and Ut- For 
T we obtain the relative error by dividing the error in T by T. For Ut we divide the error by 
the maximum of | C/" 1 1 along the curve. For both of these relative errors we take the average over 
the forty SAWs. 

We begin with the effect of the number of points we use along the curve. We study this 
using the vertical slit map with the power series approximation. We let N = 10 6 /n denote 
the number of points on the SAW used in a particular computation. For the vertical slit map 
with the power series approximation, N ranges from 10 4 to 10 6 . In figure [3J we plot the relative 
error in T as a function of N. This is the higher of the two sets of points. In figure H] we plot 
the relative error in Ut as a function of N. Again, this is the higher of the two sets of points. 
In both of these plots the error is computed by treating the values for iV = 10 6 as the exact 
values. These plots give an idea of the size of the error when the number of points used on the 
curve is reduced. 

To study the effect of the choice of conformal map (vertical slit vs. tilted slit), we compute 
the difference between the values of T using the two different maps with the same value of N. 
We do the same for Ut- We convert these errors to relative errors by dividing by T or the 
maximum of \U t \- These relative errors are plotted as functions of N in figures [3] and HI In 
both figures they are the lower sets of points. The figure shows they go to zero as iV — > oo. 
Note that these points lie well below the points that show the difference between the vertical 
map computation at the given N with the vertical map computation with N = 10 6 . In other 
words, the effect of the choice of the conformal map is smaller than the effect of using more 
points along the curve. 

Finally we discuss the effect of the power series approximation. We use the vertical slit 
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Figure 3: The top set of points is the relative error in T as a function of the number 
of points iV used on the SAW. The bottom set of points is the relative difference 
between T computed using the tilted slit and vertical slit maps. 



N 


T 


U T 


10,000 


3.66e-09 


9.12e-09 


20,000 


2.59e-09 


7.98e-09 


50,000 


1.61e-09 


4.73e-09 


100,000 


8.38e-10 


3.15e-09 


200,000 


6.98e-10 


1.76e-09 



Table 6: The effect of the power series approximation. The table shows the relative 
differences in T and in Ut when we use and do not use the power series approximation. 



method. We compute the differences in T and Ut computed with and without the power series 
approximation. We convert these differences to relative differences by dividing by T or the 
maximum of \Ut\- Table [6] shows the relative differences. They are very small and insignificant 
compared to the differences that we see when we change the conformal map used or the number 
of points along the curve. In this study we took the number of terms in our power series to 
be 12 and the parameter L = 4.0. Increasing either of these improves the accuracy of the 
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Figure 4: The top set of points is the relative error in Ut as a function of the number 
of points N used on the SAW. The bottom set of points is the relative difference 
between Ut computed using the tilted slit and vertical slit maps. 



approximation while slowing down the computation. We used block lengths of 100. 



5 Conclusions 

In this paper we computed the stochastic driving process of several models of random curves 
which we know are not SLE. We considered several statistical tests of whether this driving 
process is a Brownian motion. Simply checking if the distribution at a fixed time is Gaussian 
was seen to be useless. We must use a test that involves the independence of the increments 
of the Brownian motion. Our most effective tests were \ 2 goodness of fit tests in which we 
consider n equal increments of the process and take the cells to be determined simply by the 
signs of the increments. This test was the most successful at concluding that for the models 
that are not SLE, the stochastic driving process is not a Brownian motion. One nice feature of 
this test is that it does not involve the value of k. 

In the models we studied which are not an SLE, we have broken conformal invariance in a 
drastic way. In particular, these models are not locally isotropic. The models from physics that 
have been recently studied as possible SLE's should be locally isotropic. So a lack of conformal 
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invariance in these physical models would have to arise in a completely different way compared 
to the models we studied. Our main conclusion is not that the specific tests of Brownian motion 
that we found effective for our models are the best tests for all models, but rather that when 
we test for SLE by computing the Loewner driving process, we should employ a variety of tests 
that this process is a Brownian motion. 

We have also studied the numerical problem of computing the driving function of a given 
curve using the zipper algorithm. We have seen that the difference in the driving function when 
we use vertical slits or tilted slits for the elementary conformal map for the algorithm is quite 
small. Given that the vertical slit map is considerably faster and easier to implement, we see no 
reason to use the tilted slit map. We have also shown that the speed of this algorithm can be 
increased dramatically using power series approximations of certain analytic functions. The loss 
of accuracy from this approximation is extremely small, in particular it is insignificant compared 
to the effect of changing the number of points used to define the curve we are unzipping or 
compared to the difference between using vertical slits or tilted slits in the algorithm. 

A Details of the simulations 

In this appendix we give some details of the simulations of the distorted LERW, SAW and 
percolation interface. 

The LERW walk that we simulate is chordal LERW in the half plane from to oo. This 
means that we take an ordinary random walk beginning at the origin and condition it to remain 
in the upper half plane. Then we erase the loops in chronological order. The ordinary random 
walk conditioned to remain in the upper half plane is easy to simulate since it is given by 
a random walk beginning at with transition probabilities that only depend on the vertical 
component of the present location of the walk. If the site has vertical component k, then the 
walk moves up with probability (k + l)/4k, down with probability (k — l)/4k, and to the right 
or left with probability 1/4. (See, for example, section 0.1 of [9]). This process is known as the 
half plane excursion. The half-plane excursion is transient, i.e., each lattice site is visited by 
the excursion a finite number of times. This implies that the loop erasure makes sense. (For a 
transient walk all parts of the walk would eventually be part of a loop and so would be erased.) 
Note, however, that if we take an infinite half plane excursion and only consider the first n 
steps and loop-erase this walk, the result does not completely agree with the loop-erasure of 
the full infinite excursion. A site which is visited by the excursion before time n may be erased 
by a loop formed after time n. 

In practice there is no way to know if a visit to a site will be erased by some future loop 
without simulating the entire excursion. So in the simulation we do the following. We generate 
a half-plane excursion, erasing the loops as they are formed. We stop when the resulting walk 
has N steps. If n is small compared to N, then the distribution of our walk for the first n steps 
is close to the true distribution of the first n steps of the LERW. We will only compute the 
driving function for the first n steps. 
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We take N = 50, 000 and generate 100, 000 samples. If we work on a unit lattice, a LERW 
with N steps has a size of order N v with v = 4/5. So we rescale our walk by a factor of 
N u to obtain a curve whose size is of order one. We then compute its driving function up to 
time T = 0.01. The time T is one half of the capacity. So the number of steps needed to 
reach T = 0.01 is random. For this choice of T, the mean of this random number of steps is 
approximately 8200, roughly a factor of six smaller than N. 

The SAW in the upper half plane is defined as follows. Let N be a positive integer. We 
consider all nearest neighbor walks with N steps in the upper half plane which begin at the 
origin and do not visit any site more than once. We put the uniform probability measure on this 
finite set of walks. We let N — > oo to get a probability measure on infinite self-avoiding walks on 
the unit lattice in the upper half plane. Then we take the lattice spacing to zero. We simulate 
the SAW in the half plane with a fixed number of steps with the pivot algorithm, a Markov 
Chain Monte Carlo method |12j . We use the fast implementation of this algorithm introduced 
in [7]. For the SAW there is an issue similar to the LERW. The pivot algorithm produces the 
uniform distribution on the set of walks with N steps. But this is not the distribution of the 
infinite SAW in the half plane restricted to walks of length N. As with the LERW, we address 
this problem by simulating walks with N steps but then computing the driving function for 
only the first n steps where n is much smaller than N. 

We simulate SAWs with 200, 000 steps. We sample the SAW from the pivot algorithm 
every 100, 000 time steps in the Markov chain. We run the chain for 10 10 iterations to produce 
100, 000 samples. Unlike the other two models, these samples are not exactly independent, but 
the large time interval between sampling makes the samples very close to independent. We 
rescale the SAW by a factor of N u with v = 3/4 and then compute its driving function up to 
time T = 0.002. The mean of the number of steps needed to reach T = 0.002 is approximately 
9350, roughly a factor of 20 smaller than the total number of steps in the SAW. 

The percolation model we study is site percolation on the triangular lattice in the upper 
half plane, but we describe it using the hexagonal lattice in the upper half plane. Each hexagon 
is colored white or black with probability 1/2. The hexagons along the negative real axis are 
white and those along the positive real axis are black. This forces an interface which starts 
with the bond through the origin between the adjacent differently colored hexagons on the real 
axis. This interface is the unique curve on the hexagonal lattice which begins at this bond and 
has all white hexagons along one side of the interface and all black ones along the other side. 

Note that unlike the LERW or SAW there is no finite length effect for percolation interfaces. 
If we generate interfaces with n steps, they have exactly the same distribution as the first n 
steps of interfaces of length N where N > n. We generate 100, 000 samples of interfaces with 
N = 40,000 steps. We rescale our walk by a factor of N v with v = 4/7 and then compute 
the driving function up to time T = 0.1. This corresponds to a mean number of steps of 
approximately 11,300. 

We end with a comment on the time T and our rescaling of the various curves. For each 
model we have rescaled the curves by a factor of N u . This is merely for convenience. We could 
have left the curves on a unit lattice and computed the driving function up to time T given 
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by the above values times N 2v . What is important is that the mean number of steps of the 
curves we are unzipping is large (so that we are close to the scaling limit) but still significantly 
smaller than the total number of steps in the curve (for the SAW and LERW) so that we avoid 
the finite length effects discussed above. We have chosen the values of T so that in all three 
models the mean number of steps unzipped is on the order of 10, 000. Almost all of the time in 
these simulations is spent on computing the driving functions. Generating the random curves 
takes essentially no time by comparison. 
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